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ABSTRACT 


This report describes the results of an effort to 
refine and validate two digital microwave landing system theore- 
tical models developed under NASA Contract No. NASl-1192. The 
MLS models being considered are generic models for the Doppler 
and scanning 'beam frequency reference versions of the MLS. These 
models represent errors resulting from both system noise and 
discrete multipath. The data used for the validation effort were 
obtained from the Teocas Instrument conventional scanning beam and 
the Hazeltine Doppler feasibility hardware versions of the MLS. 
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1.0 


Summary and Suggested Studies 

A program is in progress to develop an advanced universal 
Microwave Landing System CMLS) which is intended to eventually re- 
place the present-day Instrument Landing System. The MLS consists 
of azimuth and elevation antenna systems, and Distance Measurement 
Equipment (DME) at the runway and corresponding airborne receiving 
equipment to provide position and velocity navigation data for 
terminal area flight operations. Insofar as Conventional Take-off 
and Landing (CTOL) operations are concerned preliminary design and 
system integration have been carried out by RTCA[;7] . This pre- 
liminary design included siting arrangements of the ground 
equipment, volume of coverage and format of the signals, and 
accuracy standards for the combined ground-based and airborne 
equipment. Although an effort has been made to specify the maxi- 
mum range and angle error, the statistical characteristics of the 
MLS errors are basically unknown. 

Computer simulation models that represent both the 
system noise and discrete multipath MLS errors were developed and 
implemented under NASA Contract NASI- 11992. Generic models were 
developed for frequency-reference scanning -beam and Dopple r systems. 
These models were developed from theoretical considerations con- 
cerning the communication channel. 
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The generic model developed for system noise was 
represented in terms of an autocorrelation function which was 
given in discrete form by 

(J>Cm) = e-aT|m| 

For the purpose of the modeling effort under NASl-11992, r was 
taken as the standard deviation of the elevation or azimuth error 
as given by the RTCA specifications for the MLS, and the corre- 
lation time (a was varied to determine the effect on aircraft 
performance. 

The model for discrete multipath errors resulted from 
applying Green's Theorem to determine the field arising from an 
electromagnetic wave which is reflected and scattered from an 
object such as a hangar or aircraft. The resulting discrete 
multipath errors are a function of the geometry of the multipath 
configuration as presented in Section 3.7. 

The objective of the present contract was to validate 
and refine these models by comparing and correlating the computer 
generated MLS signals with MLS signals measured from flight data. 

The data used in the study were obtained on the Texas Instrument 
conventional scanning-beam and the Hazeltine Doppler systems. 

Figure 1 illustrates the test site locations for the Hazeltine and 
Texas Instrument microwave landing systems at Wallops Island. The 
test flight path from which the data was obtained was for a straight 
in approach with a 3° glideslope. 

For the purpose of this validation effort, MLS errors 
were considered to consist of system noise, which could be con- 
sidered as a stationary random process, and discrete multipath 
errors, which result in a nonstationary random process. The 
experimental data for the system noise were used to obtain parameter 
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Figure 1. Wallops Station MLS Test Site Locations. 




estimates and these were compared with the theoretical model. 

The best fit models obtained from analyzing the experimental data 
resulted in the same form for the noise error model as that ob- 
tained from theoretical considerations, and this is considered 
to provide sufficient validation for the form of the generic error 
model which is 

♦ Cm) - 

Minor adjustments were required in the values for F in order to 
provide agreement between the theoretical models and those re- 
sulting from analyzing the data. In addition, the data analysis 
resulted in values for a which were not previously available. 

These results are discussed in Section 3.6. However, it should 
be cautioned that the results obtained on this effort were limited 
by the data that was available at the time. As more data is made 
available, the techniques presented in this report should be used 
to further refine the models. 

In the case of discrete multipath, data were not avail- 
able on the Hazeltine or Texas Instrument types to obtain a time 
history of the error characteristics. The only time history plots 
which could be located for discrete multipath errors were compared 
with the theory and similar results were obtained. These com- 
parisons are discussed in Section 3.7. Time history plots were 
used to evaluate the discrete multipath error models and favorable 
comparisons were obtained. 

As a result of the work performed on this contract and 
other recent developments in the MLS program, it is recommended 
that the following additional study efforts be performed: 

• Mathematical models should be developed and validated 
for the time^reference- scanning beam system which was 
selected by the MLS committee. 
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• The developed models should be tested with conventional 
aircraft landing systems to determined the sensitivity 
of the autopilot to correlated noise. In addition, any 
proposed autopilots developed to take advantage of the 
MLS accuracy should also be tested. 

• Optimal onboard detectors can how be developed based 
upon the spectrum of the noise statistics. Previously, 
designers could only make assumptions about the charac- 
teristics of the noise. However, since more information 
is now available, realistic designs are possible. 

• Error statistics, generated from additional airports, 
should be subjected to the same analysis. This will 
determine how sensitive the parameters of the error 
model are to site variables. At this point, all the 
software necessary to perform such an analysis has been 
developed and additional studies can be made quite easily. 
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2 . 0 Appro-ach 

The MLS system noise can justifiably be considered a 
stationary random process, and can best be analyzed in a statisti- 
cal fashion. Ultimately, an autoregressive moving- average model 
is desired. It has been shown that a model in that form can 
represent a random process using less components than either an 
autoregressive or moving- average process. However, estimating 
the parameters of such a model is complicated. The major effort 
expended in this contract was directed toward arriving at a best 
estimate for these parameters based upon the experimental data 
and comparing them with the assumed parameters of the implemented 
model. The parameters were thereby verified by using the experi- 
mental data. 

In the case of a nonstationary random process, as in 
discrete multipath reflections, a statistical verification of the 
implemented models via autoregressive moving-average models is 
ill-advised. The time histories of the error can be compared, 
however, to determine if the implemented model is reasonable. 

The procedure that was followed is outlined below; 

Cl) Develop the software to read the data tapes supplied 
by Texas Instrument and Hazeltine. 

(2) Plot the time histories of the errors in range, 
azimuth and elevation to determine any obvious problems with the data. 

C3) Develop the software to find the autocorrelation 
function and smoothed power spectrum. Plot the results for various 
values of lag, and determine the optimal lag necessary for analysis. 

C4) Develop the software to curve- fit the smoothed power 
spectrum data with a meromorphic function. 

C5) Determine the best model for the power spectrum and 
thus specify jHCj w).|^ . 

C6) From the estimated |HCjm)|^ determine the difference 
equation to model the system noise. 
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The remaining sections o£ this report will elaborate upon each of 
the six areas and present the results of the analysis. 
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3.0 


Discussion 


This section will elaborate upon each o£ the six areas 
and present the results o£ the analysis. 

3.1 Tape Read So£tware 

The £irst major e££ort per£orraed on the contract was 
directed toward developing computer programs to read the merged 
magnetic tapes which contained measured data on the Hazeltine and 
Texas Instruments £easibility hardware versions o£ the MLS and 
radar tracking data. These programs as well as the statistical 
programs were developed on the LRC Kronas Operation System using 
standard system subroutines. The tape read programs read the tapes 
and converted the BCD data to the correct £loating point value £or 
printing and plotting. The tape read programs are listed in 
Appendices A and B. 

3.2 Time History Plots 

Time histories o£ the measured errors associated with 
both the Hazeltine and Texas Instruments MLS were plotted and 
examined to determine whether there were any obvious problems with 
the data. Resulting plots o£ these time histories £or the Hazeltine 
and Texas Instrument systems are shown in Figures 2 and 3, respec- 
tively. Re£erring to Figure 2, it is apparent that £or the 
Hazeltine system, both the elevation and azimuth angle data have 
relatively large spikes which could indicate a system mal£unction 
or other large error sources in the environment. The elevation 
error spikes are particularly signi£icant because of their large 
amplitudes. It was £elt that the large error spikes should not 
be attributed to system noise; and thus, it was considered appro- 
priate to remove them Ct>y limiting the maximum error) be£ore the 
data were analyzed. The altered time history £or the elevation data 
is shown in Figure 4. Although the spikes shown in Figure 4 are 
considerably smaller than those shown in Figure 2, they are still 
present and could a££ect the analysis results. 
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3.3 Computation of Power-Spectral Density 

The next step is to compute the power spectral density. 
The spectral density of the function X(t), where X(t) is the MLS 
data, may be defined in the following manner. 

Let X^(t) = X(t) for -T < t < T 

and X,pCt) = 0 for all other values of t 

Let X^(jco) « the Fourier transform of the function X,p(t), 
which is defined as follows. 


X^Cjw} 



CtDe"^“^dt 


- CO 


CD 


■/ 


X(t)e"^“^dt 


The spectral density S(w) of the function X^t) is: 


. _ lim 1 


Xj. C-j o)) 


( 2 ) 


(3) 


If the signal XCt) is not considered. a deterministic function, but 
merely, one member of the ensemble which comprises a random process, 
the concept of power spectra must be employed to provide a harmonic 
representation of the function. Strictly speaking, such a repre- 
sentation is valid only when the random process may be said to be 
both stationary and ergodic. Briefly, this means that the statis- 
tics of the process are independent of time (no change in the 
mechanism of generation is present) and each sample function is 
representative of the whole ensemble. 


It should be noted that if the random process is not 
stationary, a power spectral representation is invalid and, in 
fact, no general representation of reasonable utility exists. If 
only the ergodicity condition is violated, the power spectrum 
representation can be used. 
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Figure 2. Time History o£ Errors for Hazeltine Data. 
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Figure 4. Corrected Time Plot of Hazeltine Data. 




A basic limitation of spectral estimation is that data 
input must always be finite in length. This causes a frequency 
smearing, or lack of resolution. This can be shown by limiting 
the data signal, X(t), by using a spectrum window function, 
as follows: 


X^Cjw) 




(t) XCt)e'^“^dt 


C4), 


where defined in several ways; the important point 

being that 

^T/2 ® ^ ^ ^MAX ^ ^ '^MAX 

= the spectral window otherwise 

From the Fourier frequency convolution theorem, -equation (4) may 
be rewritten as 


X^(w) = 


oo 

f “t/2 

- oo 


Co) - e) de 


( 6 ) 


where FC6) is the true transform and the transform of U^y 2 Ct) is 
given by 




_T . sln(MT/2') 
2tt uT/2 


if U,py 2 Ct) is the boxcar data window. Equation C6) is simply the 
convolution of two functions in the frequency domain which is 
equivalent to multiplication of the transformed functions in the 
time domain Cequation 5). 

From this it can be seen that X,p(w) is the weighted 
average of the values of F(w) about a)=6. Also X^Cu) is an estimate 
of X(w), the true Fourier transform of XCt). Because of the duality 
of time domain, multiplication and frequency domain convolution, the 
finite transform X,pCw) at u=6 is an infinite sum of contributions 
selected from X(o3) by The magnitude of these contributions 

is dependent upon the lobes of U,p ^2 either side of the maximum. . 
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It is thus desirable to minimize the size of the sidelobes of 

ill order that X^(oa) may approximate XCw). For larger data 
samples, (T-^”) , the approximation improves. In Figure 5, the 
spectrum windows used in this analysis are shown. The program 
used to calculate the spectrum generates estimates based upon both 
windows. All analyses that follow, however, are based upon the 
smooth spectrum using a Hanning window. The program used to cal- 
culate the spectrum estimate is called ASA and is available in 
the NASA Langley Scientific Library. 

Various correlation times or lags in the autocorrelation 
function were evaluated to determine the number of lags to use in 
generating the spectral density of MLS error. Figure 6 shows the 
smoothed spectrum for the elevation data, and Figure 7 the azimuth 
data, from the Hazeltine tapes. From this it can be seen that 
choosing too small a lag tends to smooth out many of the variations. 
For this reason it was decided to use 100 lags for all data. As 
can be seen with the range and angle data, very little is gained 
in increasing the lags above 100. At 150 lags the range power 
spectrum is only slightly more detailed than at 100 lags. 

Figures 8 and 9 shows the power spectrum for the Texas 
Instrument elevation and azimuth data. In this case, increasing 
the number of lags also has little effect upon the spectrum. It 
should be emphasized that the Texas Instrument tape had far fewer 
data points than the Hazeltine tapes. The confidence that one can 
have in the estimated spectrum decreases as the number of data 
points decreases. It is not wise, therefore, to attempt to achieve 
great accuracy with large lags because one is greatly limited by 
the small amount of data available. 
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Figure 6, Power Spectrum Density Hazeltine Elevation Error Lag Parameter Plots 
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Figure 7. Power Spectrum Density, Hazeltine Azimuth Error Lag Parameter Plots 












5.4 Smoot lied Power Spectra 

The next step in the analysis is to approximate the 
spectral density function by meroraorphic functions. There are two 
major reasons for using meromorphic functions. 

(1) The mean- squared error can be calculated if the 

spectral density of the error is given in the form 
of a meromorphic function of u. 

C2) Many methods exist for determining a transfer 

function which yields a minimum mean-squared error; 
these methods also are based on the assumption that 
the spectral density is a> meromorphic function of w. 

Therefore, for the application of both of the above mentioned 
methods, it is necessary first to approximate the spectral density 
curves by meromorphic functions. In addition to the above appli- 
cations, the approximation of a given curve by a meromorphic func- 
tion is of value in solving many problems involved in the synthe- 
sis of linear dynamic systems. Therefore, the following methods 
that are presented have wide application. Since this aspect of 
the analysis is in some respects the most difficult and at the 
same time the most important, a brief theoretical background dis- 
cussion is presented. 

Let us assume that we have a given curve F Cw) and that 
we wish to approximate it by a rational fractional function SCw) 
in the form 

b , + b t • • • + b 

SCu) = -2 i— 2-^ C7) 

a + a to- 1 • • • + a_to^ 

0 1 

We can write 

0 = Cb + b to^ • • • + b to^^) * S’Cto} (,3- + a +•♦•+ a ( 0 ^^) ‘ 

By selecting a? number of arbitrary points on the given curve for 
Cn+m+2) particular values of to and determining the ordinates of 
the curve, i.e., the values of the function SCto) at these points. 
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we obtain a system o£ Cn+m+2) simultaneous linear equations for 

n+m+2 unknown coefficients b , bi’*‘b„; a , a +*‘*a having 

0 a m’ ’ n ® 

the form 


0 


(I’o * I’l"/ 


SCw^) 



+ a w ^ 
a 1 


+ • 


• +a„w 
n j 



^ ^“n+n,*2 

a 

n n+m+ 2 y 


C9) 


CIO) 


The solution of these equations will provide us with the numerical 

values of the coefficients b , b„b ^ ^ . 

0 m m a n 

Besides being awkward, this method has the further dis- 
advantage that even if the function SCw) is continuous and the 
number of points taken is very large, it is still impossible to 
verify that the error C^) in the approximation, i.e., 

A = F(w) - SCw) (11) 


is of sufficiently low absolute value at points which do not 
coincide with the selected points u) 2 »**. 


An alternate approach to solving this problem is using 
nonlinear least squares. The standard method for solving least 
squares problems which lead to nonlinear normal equations depends 
upon a reduction of the residuals to linear form by first-order 
Taylor approximations taken about an initial or trial solution for 
the parameters. The nonlinear nature of the problem can best be 
seen if SCw) is written in the following form: 


S(w;a,b) 


a + a + • • • + 
0 -1 


b . + b t • ' » + 
0 . 


n-j 

b 0)2^ 
n 


( 12 ) 


where the a’s and b's denote the numerator and denominator coefficients, 
respectively. The least squares solution to this problem by 
direct minimization of the approximation error as a function of 
the filter coefficients requires the solution of a set of non- 
linear equations. 
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Several methods 11-3] are available which have decoupled 
the solution for the numerator and denominator coefficients. First 
a linear estimate for the b*s is calculated and then the corres- 
ponding a’s are found. However, these linear estimates are sub- 
optimal, since the b’s are solutions to an overdetermined set of 
linear equations that minimize the linear equation error eCh). 

This e(b) is a nonlinear mapping of the approximation error e(a,b). 
In general, the solution for b that minimizes -| IeCb)l| is only the 
same as that which minimizes ||eCa,b]|| when eCa-,b) = 0. Therefore 
it is necessary to minimize with respect to both a and b. 

A program has been used which solves for a and b simu- 
taneously. The method is based on the Marquardt technique 14] 
which is similar to the Levenberg procedure £5]. 

^in - 

Let (j>CB] -^[W. Y. - FCX^,B)]^ C.13] 

i=. i 

The usual linearized Gauss procedure [6] applied to the problem 
of minimizing i()(B) generates estimates/ by the following procedure: 

Start with an initial estimate B. 

Solve . 

a"^A 6 = A*^r for 6. C14] 

The new estimate is calculated as B^ = B + tS where 0 < t < 1. The 
process is then repeated at the new estimate B^ until some accep- 
table convergence criterion is met. In the above, 

r. = the ith residual = W- - FCX^,B)^ 

A = the matrix defined by 

Aji^j = I partial of FCY-,B) with respect to bjj 

6 = the correction vector, B^ = B + td. 

T T 

The Marguardt method differs from the above in that CA A+^) 6= A x 

T T 

is solved instead of A A = A r . The J value is called a damping 
coefficient. In practice, the method is implemented in the pro- 
cedure given below. 
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Start each iteration with a Bo and a ^o.. Bo is the 
current estimate of the minimizing parameter vector and Xo is an 
estimate of the "scaled” damping coefficient. 

X T 

The A A and A r quantities are formed at the value Bo. 

This system is scaled to form {Sa'^AS + XoS^jS~^6 = SA^r. Thus, the 

Xo above is an estimate of the damping constant for this scaled. 

T -1/2 

system. Here S is a diagonal matrix with S. ^ = lA A. -3 ' . The 

scaled correction factor, S "^5, is computed and the new correction 
is Bi-= Bo + ^ 6) . -If- f (Bj) < then a new vector is calcu-^ 

lated using X/io, this is called B 2 . The new X estimate is chosen 
to be the one corresponding to minimize ^CB 2 )J • The cor- 

responding Bj^;K=l,2 is taken as the next parameter estimate. If 
<j>CBj^) > (j>CBo), the new vector, B 2 , corresponding to a damping con- 
stant of X/10 is still calculated. If 4>CB2D <}>CBo)» X/10 is 
accepted. If (j)CB 2 ) > 4>CB2)» then X is increased by factors of 10 
until one of two conditions occurs. If some X*10^ is found that 
yields an estimate that decreased <J), this X* 10 is taken as the 
damping constant and the new estimate is the B 2 found from 
CSA^AS + X*10^S^) S~^6 = SA'^r. However, if the angle between the 
Marguardt direction vector and the negative gradient direction 
becomes smaller than 45 degrees before such a X*10 is found, then 
the search on X is terminated and the direction taken to be examined 
is the negative gradient. 

The entire procedure is repeated until the maximum 
number of iterations is exceeded or convergence is achieved. 
Convergence is accepted when -|6. | < r+e|B^| for i*l,no. Here 6 
is the correction vector to B and r<<e. 


Since the Marquardt method requires many calculations, a 
linear least squares method was developed and compared to the 
Marquardt method. 


Let 

SC“) 


aa +. a + a w 

■ 0 X n 


1 + b • • • + b^^ 

. a n+a 


zn 

InTTT 


CIS) 


23 



This can be written as 

S (to) + b S(to)to^ +••*+ SCw)b ,(0^^^^^^ a + a (o^ +»••+ a to^’^ 
4 n +1 a 4 n 

S (to) = S(to)b (0^ - SQo)b , + a + a to^ +•••+ a to^^ 

4 n+± ® 1 n 

Making a change o£ variables as follows: 



and 

= -SCto)to^ 

Y2 = -SCto)(o^ 

• * 

Y^ = -S(to)(0^'^ 

Results in 

SCto) = b Y + b Y +•••+ b„Y„ + a + aX + aX +••«. ( 16 ) 

The transformed equation is now linear in its coefficient wLich. 
implies that the "a's" and ’*b’s” can be estimated by linear regres- 
sion. The nonlinear correlation between to and S(to) will not 
affect the efficiency of estimation. In this form, the estimation 
problem has been greatly simplified. A computer study was made 
to verify the performance of the linear regression as compared to 
the nonlinear regression. In general, th.e linear regression tended 
to give estimates which weighted all data points equally; whereas. 


24 



the nonlinear program tended to weigh more heavily those data 
points of high magnitude. It is felt that additional work should 
be done with the linear regression because theoretically it repre- 
sents a great simplification of the approximation problem. Indeed, 
it is possible to estimate the coefficients on-line so that one 
could get a fast estimate of how quickly the error statistics are 
changing. This method could also be used in designing digital 
filters where the impulse response is given and the transform HCZ) 
is desired which will best approximate the impulse response. 

3.5 Best Fit Models 

The math models which give the best fit to the Spectra, 
based upon the nonlinear regression program are given below; and 
are illustrated in Figures 10 through 14. 

Hazeltine Data: 

C17) 

0.4288 
28.02 

C18} 

0.01314 
0.128 

C19). 


Elevation 


HCjw) 


B 


CB +w^) 
1 


B = 7.024 Standard Error = 
0 

B = 366.37 Standard Error = 
-1 


Azimuth 

iHCjw) 


B 


CB +03 } 
1 


B = 0.1156 Standard Error = 
0 . 

B = 0.938 Standard Error = 
a 


Range 




B 


Cb~ho 


B = 1.4998 Standard Error = 2.62 

B = 1.02389 Standard Error = 0.1092 
1 

Texas Ins trumeht : 


B 

0 


Azimuth 

|HCjo3) 1^ 


C20) 



lHi(ja)>|**.Q5 (DEG) 
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0.004 
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Figure 10. Hazeltine 
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Figure 11. Hazeltine Azimuth Data. 
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Figure 12. Hazeltine Range, Power Spectrum Density. 
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Figure 13. Texas Instrument Elevation Data. 
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Figure 14. Texas Instrument Azimuth Data. 
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B = 0.33025 Standard Error = 0.026 
= 0.567 Standard Error = 0.051 
Elevation g 

V (21) 

CB +0) ) 

1 

B = 0.687 Standard Error = 0.0634 
0 

B = 5.066 Standard Error = 0.5666 

X 

Several additional models were fit to the data but the results in 
general were poorer. Table 1 shows a listing of the models with 
the standard errors of the coefficients. With several models, 
the coefficients did not converge to a final value; this indicates 
that the hypothesized model is an inadequate representation of the 
data. It should also be pointed out that the models that did the 
best for angle data, for both the Hazeltine and Texas Instrument 
data, are also the simplest with which to work Cho zeros and just a 
single pole). • The specification of the difference equation becomes 
quite straightforward and the results are more easily understood. 
For instance, it is possible for the angle data to specify the 
errors by the correlation time, a variable that is derived simply 
from the "B” coefficients. This simplicity is a very desirable 
property, one which cannot be too strongly emphasized. 

An important fact to observe from the time plots of the 
Texas Instrument data (Figure 3) is the small time period of the 
available data. The plot of the azimuth angle error data indicates 
that its statistics are not stationary. This is also true for the 
range data where the first. four or five seconds of data are ob- 
viously nonstationary. Therefore, a large amount of confidence 
cannot be placed on the power spectral estimates for the Texas 
Instrument data. As mentioned earlier, power spectral estimates 
are not valid for a nonstationary process. In addition, it was 
extremely difficult finding transfer function coefficients that 
would represent the data and the number of data points available 
for angle analysis both for azimuth and elevation data was less 
than 200. This number is far too small for adequate analysis. 

For these reasons, it was felt that to carry the Texas Instrument 
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angle data analysis any further would be misleading. Basically, 
the data does not lend itself to the type of analysis being 
performed, and to continue with the analysis would be misleading. 

The range data, excluding the first five seconds, are 
also open to question. Unlike the Hazeltine data, there is clearly 
a bias on the range error data (compare Figures 2> and 3^). An 
attempt was made to continue the analysis with the first five 
seconds removed from consideration. Even then, however, it must 
be pointed out that the results obtained may be misleading. 


5'f & i V a tionPi iff e'/r erit ®qiuat i-ons. c£o>ri sDigital Simulation 

A detailed derivation of the Hazeltine elevation 
difference equation is given. Since the derivation for the 
other models is similar, only the results are presented for them. 

Hazeltine Elevation: 



9 ® 0 

Bq = 7.024 


H(jw) 

1 ■ ■ 

' 

C22) 

'.(B + u^) 

1 

Bi = 366.37 


H(jo)) 

2 _ — 1y-2u- . 
' 0 . 2 -^*- 0)2 

a = 19.1 
Y = 0-.a84 

(23) 


^(t) = Autocorrelation function 
= Scale Factor e •' 


The scale factor results from the normalization that the sub- 
routine ASA performs. In general = mean-squared error 





EIX(f)^J 


[ECX)J.- + Variance 


(24) 


Therefore, in order to properly scale the autocovariance, E(X) and 
2 

o must be known. It came to our attention at the conclusion of 
this project that the subroutine ’'ASA*' which calculates the power 
spectrum density calculates the following parameter. 


2 , X- E jjC] 
a 
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and calculates the spectrum for Z and not X. This in itself is of 
no major concern as long as ECX) and a are both kno;vn. However, 
a is not printed by ASA; in fact, the entire normalization pro- 
cedure is never mentioned in the documentation. 


The normalization procedure does not affect the shape of 
the spectrum, but it does affect the scale factor. Therefore, 
all results for the mean-squared error are conservative (i.e., the 
results do not include the variance which is unknown). Continuing 
in light of the above, the autocorrelation function is given as 
follows: 

$(t) = C.184) 

$(t) = 50.703 X 

= r^e~“ I 


r = 7.01 X 10'^ deg 
sec ^ 


In discrete form, 

tCm) - r2e-“Tl”>r 


(25) 


The 1 transform of $Cni) is obtained by taking the sum of the 
individual Z transforms of the parts for m > 0 and m < 0. Letting 
A=e-“T 


then 


we have 


$CZ) 

-{ ^ .1 ^ 

*(Z) 

■f'^1 


\l-AZ > J 

H(Z) 

ri/r-A* 


1-AZ‘ ^ 

H(Z) 

= : [1 

fi-CO). 

= rc/r^l 

^00 

: /$t0) (01 




C26) 

C27) 

(28) 

(29) 

C30) 

( 31 ) 
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Cq = CqqVCO) = TAVCO) C32) 

yCO) = r/l-A^ uCO) + TAVCO) (33) 

Since u(0) and v(0) are independent and tKeir values do not enter 
the expression for yCn) for n, > 1, y(0) can be generated more simply 
from a single random variable having the appropriate variance, by 
taking 

y(0) = uCO) C34) 

Finally for n > 1 

yCn) = r/l-A* liCn) + AyCn-1) C35) 

Since F represents the mean- squared error, this variable 
can be compared to the total error that is specified by RIGA. A 
copy of this is reproduced in Table 2, As can be seen the 
Hazeltine system can be used for configuration k. Cat in, based 
on azimuth data. The remaining results for the Hazeltine data 
are shown in Table 3. 

Several significant conclusions can be made about the 
results presented in Table 3. First, the correlation time (i) is 
considerably different for azimuth and elevation data, varying by 
a factor of twenty to one. Also, the correlation time for range is 
similar to the azimuth. The mean-squared error for range is just 
barely up to the requirements, and in the case of elevation, the 
requirement has not been met. 
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Table 2. 

RTCA MLS, Specifications. 
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Random error negligible compared to bias. 



Hazeltine MLS Error Models. 




3.7 


Discrete Multipath Errors 

The preceding discussion has described the validation 
results for the system noise component of the MLS error. The 
other error component included in the MLS model results from 
discrete multipath reflections. Data were not available on the 
Hazeltine or Texas Instrument tapes to obtain a time history of 
the multipath characteristics; however, Reference 8 contains the 
results of a study of discrete multipath effects. Figure 15 
shows this effect for the elevation beam of a time reference MLS. 
Figure 16 displays the error for the localizer of a frequency 
reference MLS based on the simulation models described in 
Reference 9. 

As can be seen for both error models the magnitude of 
the error builds up over a relatively long period of time, reaches 
a maximum and then decays. This effect can be explained by the 
model used for the discrete multipath reflections. As the air- 
craft approaches the specular reflection point the error due to 
reflections increases. Outside the specular reflection point the 
error is less but somewhat periodic. It would be highly desirable 
to develop simulation models for discrete multipath errors for the 
time reference MLS and to obtain empirical data to better verify 
this type of performance on both channels. 
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Figure 15. Scanning-Beam Angle Error, 40-Hz Sampling Rate (output from Filter). 
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Figure 16, Time Plot of Discrete Multipath Errors. 
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Figure 17. Airport Geometry. 
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IF(tFG.EQ.JFG(l)lGO TO 10 

if(abs(difr).gt.abs(pt))go to 10 
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Appendix A. Texas Instrument Fortran Extended Read Program for Data, CContinued) 
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APPENDIX C 

NONLINEAR LEAST SQUARES PROGRAM 


Note: This program was adapted from' 

"The Computing Technology Center 
Numerical Analysis Library" CTC-39 
G.W. West ley, Oct. 1970. 
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